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ABSTRACT 

Those galaxy clusters which do not belong to the superclusters are referred 
to as the isolated clusters. Their relative abundance at a given epoch may be a 
powerful constraint of the dark energy equation of state since it depends strongly 
on how fast the structures grow on the largest scale in the Universe. We note that 
the mass function of the isolated clusters can be separately evaluated through the 
modification of the recently developed Corasaniti-Achitouv (CA) theory accord- 
ing to which the stochastic collapse barrier is quantified by two coefficients: the 
drifting average coefficient (/?) and the diffusion coefficient {Db)- Regarding (3 in 
the CA formalism as an adjustable parameter and assuming that the formation 
of isolated clusters corresponds to the case of Db = 0, we determine the mass 
function of the isolated clusters by fitting the numerical results from the MICE 
simulations to the modified CA formula. It is found that the best-fit value of 
changes with redshift and that the CA mass function with Db = agrees very 
well with the numerical results at various redshifts. Defining the relative abun- 
dance of the isolated clusters, £/, as the ratio of the cumulative mass function 
of the isolated clusters to that of the non-isolated clusters at a given epoch, we 
finally show how sensitively £/ changes with the dark energy equation of state. 
It is also discussed how can help to break the degeneracy between the dark 
energy equation of state and the other key cosmological parameters. 

Subject headings: cosmology:theory — large-scale structure of universe 

1. INTRODUCTION 



The observed rich clusters of galaxies are believed to have form ed through the grav i- 
tational collapse of the highest local peaks of the initial density field ( IBardeen et al.l 119 86). 
Although most of the rich clusters belong to the larger-scale superclusters (e.g., IWray et al. 
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20061 ; lEinasto et al.ll2007| ). there exist few isolated massive clusters which reside in relatively 
low-density environments without having any close neighbor clusters. The relative abun- 
dance of the isolated clusters should depend on how strong the gravitational clustering is on 
the largest scale and how frequently the clusters merge into the superclusters in the Universe. 
Therefore, it is expected that the relative abundance of the isolated massive clusters may be 
useful to distinguish between different candidate dark energy models. 

To predict the relative abundance of isolated clusters in different dark energy models, 
however, a theoretical framework is first required within which the number density of iso- 
lated clusters can be separately counted. The standard excursion set formalism has been 
conventionally employed to analytically calculate the number density of galaxy clusters as 
a function of mass, i.e., the mass function of galaxy clusters. Basically, the standard excur- 
sion set formalism counts the initial regions whose linearly extrapolated density contrast (5) 
exceeds a fixed threshold value (S c ) on a certain mass scale M, regardin g them as the sites 
which would eventually collapse to form the bound objects of mass M (IPress fc Schechter 
1974 ; iBond et al.lll99ll ). The density threshold S c which is often called the collapse barrier 



has a cons tant value of = 1 .686 if the gravitational collapse process follows the spherical 



dynamics ( Gunn fc Gott 19721 ). 

Maggiore fc Riotto have for the first time introduced the concept of the 



stochastic collapse barrier to gener a lize t h e excursion set formalism. In the light of the 
numerical result of Robertson et al.l ( 2009 ). Maggiore &; Riotto ( 2010a b) regarded the col- 
lapse barrier, 8 C , as a log-normal variable, and showed that in the high -mass limit the collapse 
barrier has a spherical average of (5 C ) = S sc = 1.686. Very recently, ICorasaniti fc Achitouv 
( 12011al Jbl) have made a further refinement of the mass function theory by incorporating the 
ellipsoidal collapse dynamics into the generalized excursion set formalism. They extended 
the generalized excursion set formalism to the case that the mean collapse barrier, (S c ), de- 
viates from the spherical average, 5 SC = 1.686 in the low-mass limit. The higher collapse 
barrier than the spherical average, 5 C > 1.686, corresponds to the case that the formation of 
the low-mass halos experience the disturbing tidal effect from the surroundings. 

We note here that the Corasaniti-Achitouv formalism should be applicable not only to 
the case that the gravitational collapse occurs non-spherically in the low-mass section but 
also to the case that the massive halos form in isolated low-density environments. Given that 
the massive cluster-size halos are likely to form in the highly overdense regions but hardly in 
the isolated low-density regions, the collapse barrier for the formation of the isolated clusters 
should be higher than the spherical average. Since the Corasaniti-Achitouv formalism allows 
the average collapse barrier to deviate from the spherical value of 1.686, their formalism may 
be capable of deriving the mass function of the isolated clusters. 
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The organization of this paper is as follows. In §2 we briefly review the Corasaniti- 
Achitouv formalism of the halo mass function. In §3 we explain how the relative abundance 
of the isolated clusters can be evaluated by modifying the Corasaniti-Achitouv formalism 
and present the results of the numerical tests. In §4 we show how sensitively the relative 
abundance of the isolated clusters depends on the dark energy equation of state as well as 
the other key cosmological parameters. §5 we summarize the results and conclude the work. 



2. THE CORASANITI-ACHITOUV FORMALISM 

The mass function of bound halos, dN / d In M, is defined as the number density of the 
bound objects whose ma sses belong to the differential mass interval [In M, In M + d In M\ 



( jPress fc Schechterl Il974j ). The classical mass function theory based on the standard ex- 



cursion set formalism relates dN/dlnM to the multiplicity function, f(cr), which gives the 



(Bond et al. 


1991; 


Jedamzik 


1995) 



dN(M,z) p diner' 



cilnM 



M dlnM 



■f[a(M,z)}. 



(i; 



Here p is the mean mass density of the Universe and the pseudo time variable cr(M, z) is 
equivalent to the rms density fluctuation of the linear density field smoothed on the mass 
scale M at redshift z. It scales with the linear growth factor, D(z), as a(M, z) = D(z)a(M) 
with a(M) = a(M, z = 0). Throughout §2 and §3, we assume a flat ACDM cosmology (with 
Q m = 0.25, tt A = 0.75, Q b = 0.045, h = 0.7, n s = 0. 95 and a s = 0.8) for which the linear 
growth factor has the following analytical expression (ILahav et al.lll99ll ): 



D(z)^-n m [n m (i + zf + n A ) 1 / 2 



Here D(z) is normalized to satisfy D(z = 0) = 1. 



dz' 



1 + z' 



[^ m (l + z') 3 + ^A] 3 / 2 ' 



(2) 



In the original Press-Schechter theory where the gravitational collapse is assumed to oc- 
cur spherically, the collapse barrier has a constant value of 5 SC = 1.686. In the Press-Schechter 



ing for the ellipsoidal collapse process (e.g., 


Bond & Mvers 


1996 


: I 


Vlonaco 


1997; 


Audit et al. 


1997; 


Lee & Shandarin 


1998 


Sheth & Tormen 


1999 


Jenkins et al. 


2001; 


Warren et al. 


2006; 


Tinker et al. 


20081 


Robertson et al. 


2009). 


the collapse barrier is described not as a con- 



stant but as a function of the mass scale M (or equivalently, a function of the rms density 
fluctuation, a). 
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Recently, iMaggiore k, Riottol ( l2010aU bl. hereafter MR10) pointed out that not only the 
constant spherical collapse barrier but also the scale-dependent ellipsoidal collapse barrier is 
incapable of describing the true nature of the complicated halo formation process and sug- 
gested that the collapse barrier S c should be treated as a stochastic variable. MR10 derived an 
analytic expression for the halo mass function by generalizing the excursion set formalism for 
the case that the collapse barrier is stochastic and the random walk process is non-Markovian. 
It was shown by MR10 that in the high-mass limit (M > 10 14 /i~ 1 M ) the average of the 
stochastic collapse barrier equals the constant spherical value, (5 C ) = 5 SC = 1.686 and its vari- 
ance is scale-dependent as (A 2 Db) = ct 2 (M)Db where Db is called the diffusion coefficient. 
The stochasticity of the collapse barrier is caused by the di sturbance from the surroun dings 
and the ambiguity in the halo-identification procedure (see IMaggiore &: Riottoll2010bl . for a 
detailed explanation). 

In the light of the MR10 work, Corasaniti fc Achitouv ( 2011al lbl. hereafter CA11) have 
extended the generalized excursion set formalism to the ellipsoidal collapse case in which the 
average of S c deviates from the spherical value, 1.686, in the low-mass section. The resulting 
CAll multiplicity function has two characteristic coefficients: the diffusion coefficient Db 
defined as in the MR10 formalism and the drifting average coefficient, 0, defined as (5 C ) = 
8 SC + (3. The non-zero value of this drifting average coefficient, 0, quantifies the degree of 
the deviation of the average of the collapse barrier from the spherical average 1.686, while 
the non-zero value of the diffusion coefficient Db quantifies the degree of the stochasticity 
of the collapse barrier S c . 

The CAll multiplicity function f(<j) was written at second order as 

f(a; D B , /?) « f {0 \a; D B ) + ffl Q {a; D B ) + ff\a\ D B ) + f$(a; D B ), (3) 



where 



ay/1 + D B V 7r 



2 _ (<5 sc +/3<t 2 ) 2 




a4 ^ 1 ^ / 
e -^r _ -T(0,x) 



a S sc | fpi (cr) + R Erfc (x) 



(4) 

(5) 
(6) 



f^D B ) = ^ 2 ^{Erfc(,) + -J-^ 



e 2ct- 



1 n \ 3a5 2 
2x + — 

2 A 



^fr(0,a;)]},(7) 



with x = (5 sc /a)y/a/2, a = 1/(1 + Db), k, = kcl, k = 0.475, and incomplete Gamma 
function r(0,x). Here, f(°\a; Db) coincides with the MR10 multiplicity function which 
has only one parameter Db- Computing the halo mass f unctions through E quations (JTJ- 
(JTJ) and comparing them with the fitting formula given by iTinker et al.l (120081 ) . CAll have 
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determined the best-fit values of the two coefficients as (3 = 0.057 and Db = 0.294. CAll 
have also shown h ow the shape of the halo mass f unction changes with the values of Db and 
/3 (see Figure 3 in ICorasaniti &: Achitouvll2011bl ). According to their results, the variation 
of the diffusion coefficient Db alters the high-end slope of f(a) while the variation of the 
drifting average coefficient f3 affects the over-all amplitude of /(cr). 

Now, we would like to confirm the vali dity of the CAll m ass function by testing it 
against the high-resolution MICE simulations ( jCrocce et al.ll2010l ) which traced the evolution 
of 2048 3 dark matter particles (each having mass of 23.42 x 10 10 h Mq) on a periodic box 
with linear size of 3 /i -1 Gpc for a ACDM cosmology. We utilize the publicly available catalog 
of the cluster halos that were identified via the standard Friends-of- Friends (FoF) algorithm. 
The catalog provides information on each halo's (comoving) position, peculiar velocity and 
FoF mass (calculated as the sum of the masses of all the dark matter particles belonging to 
each halo) at three different redshifts: z = 0, 0.5, 1. For a d etailed description o f the MICE 
simulations and the cluster catalogs, we refer the readers to ICrocce et al.l (120101 ) . 



Binning the logarithmic masses of the cluster halos from the MICE catalog and counting 
the number densities of the cluster halos belonging to each logarithmic mass-bin, we obtain 
the numerical mass function of all cluster halos, dNx/dhiM, at each redshift. Fitting 
the CAll formula (Eqs [I]- [7]) to the numerical results from the MICE simulations at each 
redshift, we have determined the best-fit values of (3 and Db with the help of the \ 2 statistics, 
the results of which are listed in Table [TJ 

Note that the diffusion coefficient, Db-, has the same best-fit value of 0.38 at all three 
redshifts while the best-fit value of the drifting average coefficient, /3, increases as the redshift 
decreases. The result indicates that at earlier epochs the gravitational collapse process is 
closer to the spherical dynamics than at present epoch. Since a bound halo is harder to form 
at earlier epochs, those high- z halos must correspo nd to higher peaks in the initial density 
field. As shown analytically by iBernardeaul (jl994J), the higher a local density peak is, the 
more spherically its gravitational collapse proceeds. 

It is also worth mentioning that our result on the best-fit value of Db = 0.38 is different 
from the original value Db = 0.27 used in the CAll formalism. This discrepancy must be 
due to the fact that we use the FoF masses available in the MICE catalog while the spherical 
over-density (SO) masses were considered in the original CAll work. 

Figured] plots the CAll mass functions with the best-fit values of (3 and Db (solid line) 
and compare them with the numerical results from the MICE simulations (square dots) at 
z — 0, 0.5 and 1 in the top-left, top-middle and top-right panels, respectively. The Jackknife 
method is employed to calculate the errors: Dividing the halos into eight subsamples (each 
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having the same number of the halos) and determining dNx/dhiM separately from each 
subsample, the errors are calculated as one a scatter of dNx/d\a M among the eight Jackknife 
resamples. Figure [T] also plots the ratio of the analytic mass function to the numerical result 
as a function of mass at z = 0, 0.5 and 1 in the bottom-left, bottom-middle and bottom-right 
panel, respectively. As can be seen, the analytical and numerical mass functions agree with 
each other quite at all three redshifts, except in the high-mass section (M > 10 15 hr x M Q ) 
where the Jackknife errors are very large. 

3. MASS FUNCTION OF THE ISOLATED CLUSTERS 

Applying the FoF algorithm with the linkage length parameter of b to the MICE cluster 
catalogs at each redshift, we first find the clusters of clusters (i.e., superclusters) which have 
more than one neighbor clusters within the FoF linkage length, bl where / is the mean cluster 
separation. The isolated clusters are then identified as those clusters which do not belong 
to any superclusters. Note that the degree of their isolation depends on the value of b: 
The larger the value of b is, the more isolated they are. Here we consider the extreme case 
in which the isolated clusters experience no disturbance from the surrounding large-scale 
structures and the difference between their FoF and SO masses is negligible. In this extreme 
case, the value of the diffusion coefficient, Db, must vanish since the non-zero value of Db 
indicates the presence of the disturbance from the surroundings and the existence of the 
difference between the FoF and SO masses. 

The ad-hoc value of the linkage length parameter for this extreme case is determined 
to be b = 0.4 by the following procedures. We first investigate how the best-fit value of Db 
changes as the value of b varies from 0.25 to 0.45 in the FoF algorithm applied to the MICE 
cluster catalogs. For each case of b, we select the isolated clusters from the MICE cluster 
catalog and obtain the numerical mass function of the isolated clusters. Then, we fit the 
numerical result to the analytic CA11 formula (Eq.[T]-[7]) to determine the best-fit values of 
D B and (3 with the help of the \ 2 statistics. Figure |5] plots the best-fit value of D B versus 
the linkage length parameter b. As can be seen, the best-fit value of Db drops to zero when 
the linkage length parameter b reaches up to 0.4. Table [2] lists the best-fit values of Db and 
(3 for the mass function of the isolated clusters when the ad-hoc value of the linkage length 
parameter is set at b = 0.4. 

Now, we write the total mass function of all cluster halos, dN^/d In M, as the sum of the 
mass function of the isolated cluster halos, dN\/dhiM and that of the non-isolated cluster 
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halos, dN m / ' d In M: 



dN T (M, z) dNi{M, z) dN m (M, z) 



dlnM dlnM dlnM y ' 

Putting Db = in the CAll formula, we express the mass function of the isolated clusters 

as 



dNj(M, z) p dlna~ x 



f[a;(3,D B = 0}. 



(9) 



dlnM M dlnM 

Here the multiplicity function of the isolated clusters which is characterized by one coefficient 
(3 is written at second order as 



f(a; D B = 0, fi) « + A%(a) + ff\a) + f$\* 



,02 



where 



/ (0) M 



— K- 



6 sc 1 „ , 

e -^__r( 0) x) 



-ps. 



4 1 2 ( ( x) + «:Erfc(x) 



/3 2 5 s 2 c fi:<^Erfc(a;) + 




(7 

5 SC V 27T 



3 5? 



f i --2.r ! f -^r(o,x) 



(10) 

(11) 

(12) 
(13) 
(14) 



Figure [3] plots the analytical mass functions of the isolated clusters (solid lines) with the 
best-fit values of (3 and compared them with the numerical results from the MICE simulations 
(square dots) at z = 0, 0.5 and 1 in the top-left, top-middle and top-right panel, respectively. 
As can be seen, the analytical mass functions of the isolated clusters agree quite well with 
the numerical results at all three redshifts. The ratios of the analytic CAll formula to the 
numerical results are plotted in the bottom panels. Although the agreements between the 
analytical and the numerical results for the case of the isolated clusters are not so excellent 
as for the case of all cluster halos (see Figured]), the ratio at each redshift is still quite close 
to unity: the discrepancy is less than 20% except in the high-mass limit where the Jackknife 
errors are very large. 

This result indicates that the CAll formula with D B = indeed works fairly well for 
the determination of the mass function of the isolated clusters. It is worth emphasizing here 
that this analytical mass function of the isolated clusters (Eq. [5]) has only one coefficient 
f3 and there is no fitting normalization constant. The amplitude of dN\/ dlnM is found to 
automatically match the numerical result when the best-fit value of [3 is put into the CA 
formula. 
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Figure H] plots the ratio of the mass function of the isolated clusters to that of the non- 
isolated clusters at z — 0, 0.5 and 1 in the left, middle and right panel, respectively. The 
errors are calculated as one a scatter of the ratio among the eight Jackknife resamples. As 
can be seen, the numerical result (square dots) agrees fairly well with the analytic prediction 
(solid line) based on the CAll formula at each redshift. We believe that it would reduce the 
existing discrepancy between the analytical and numerical results even more to include the 
higher order terms in the calculation of the multiplicity function, f(cr; Db, 0)- Throughout 
this paper, however, we consider only the second order approximation of f(a; Db, 0). Note 
that the ratio drops with mass more rapidly at higher redshifts, indicating that at high 
redshifts most of the very massive clusters have stronger tendency to form in the highly 
overdense regions than at present epoch. 



4. DEPENDENCE ON THE DARK ENERGY EQUATION OF STATE 



In a ACDM cosmology the dark energy equation of state (defined as the ratio of the 
pressure density to the energy density) is a perfect constant given as w = -Pa/Pa = —1. 
Whereas in dynamic dark energy 



models the value of w may vary with time and can deviate from — 1 at z = ( IWang fc Steinhardt 



19981 ). This dif ference in the dark energy equation of state results in the different funct ional 
shape oi D(z) Jchevallier fc PolarskilEoOll iLinderihoOsI : lBasilako]l2003t Ipercivalll2005l ). 



The evolution of the abun dance of galaxy clusters ha s been regarded as one of the most 



powe rful probes of w(z) (e.g., IWang &: Steinhardtl Il998l : lHaiman et al.l l2001t IWeller et al. 



20021 ). The most serious systematics in constraining w(z) with the evolution of the abun- 
dance of the galaxy clusters may come from the mass assignment of the high-z clusters. 
Several statistical methods have been suggested so far to overcome the systematics but none 



of them have yet to be fully satisfactory (e.g., iRykoff et al.l 120081 ; ICunha fc Evrardl 12010 
Stanek et al.ll2010l . and references therein). See also Allen et al.(2011) for the latest review. 
We suggest here that the relative abundance of the galaxy clusters at a given epoch provides 
a complimentary probe of the dark energy equation of state. 

The relative abundance of the isolated clusters is defined as the ratio of the cumulative 
mass function of the isolated clusters to that of the non-isolated clusters: 



£i(M c 



dlnM 



dNi(M, z) 
dlnM 



d\nM 



M c 



dN m (M,z) 
d\nM 



(15) 



where the mass function of the non-isolated cluster dN-^i/ dhiM is obtained by dN^/dlnM — 
dNi/dlnM. To demonstrate how the relative abundance of the isolated clusters £i changes 
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with the dark energy equation state w, we consider the specific case where the dark energy 
equation of state i s r edshift de p enden t as w(z) = w n + W\ z/{1 + z) ( IChevallier fc Polarski 



200 ll ; iLinderl 120031 ). iBasilakosi ( 120031 ) and iPercivall (120051 ) showed that for these dynamic 



dark energy models the linear growth factor can be approximated as: 



D(z) 



5a 



2(z + 1) 



Q a mz -Q Qz+ [l + 



(l + AQ Qz ) 



where 



a 



A 



+ 



(1- W )(1-3 W /2) M , 



5-2/(l-w) 125 (l-Qw/5) 3 
0.28 

0.3. 

w + 0.08 



(16) 

(17) 
(18) 



Here Q mz and Qq z represent the matter density and dark energy density parameters at z, 
respectively, related to their present values, Q m and Qq as 



a 



Q m (l + z) 
E 2 (z) 



n 



Qz 



E 2 {z){l + z)fW 



where 



E\z) 



n m (i + z) 3 + n Q (i + z)-^ z \ 

3wi 



-3(1 + w ) 



21n(l + z)' 



(19) 

(20) 
(21) 



We first calculate £i for the ACDM case and then repeat the calculation of £i for the four 
different dynamic dark energy models at z — 0.5. Figure |5] plots the relative abundances of 
the isolated clusters for the dynamic dark energy models and compare them with the ACDM 
case. The errors for the ACDM case are calculated as one a scatter of £i among the eight 
Jackknife resamples from the MICE cluster catalog. As can be seen, the five different dark 
energy models predict different relative abundances of the isolated clusters. For the case of 
W\ > (wi < 0), the relative abundance of the isolated clusters has lower (higher) amplitude 
than for the ACDM case. This can be explained by the following logic. In a dark energy 
model with w\ > (wi < 0) the largest-scale powers are higher (lower) than in the ACDM 
case and thus boost the merging of the isolated clusters into the superclusters, reducing the 
relative abundance of the isolated clusters. 

Furthermore, the differences in £i between each dynamic dark energy case and the 
ACDM case are larger than the statistical errors calculated through the Jackknife resampling. 
If a future cluster survey can find as many clusters as the MICE simulations, the sensitivity 
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shown in Figure [5] suggests that the relative abundance of the isolated clusters will be a 
powerful probe of the dark energy equation of state in practice. 

We examine the degeneracy between wq and W\. Varying the value of Wq from —1 to 
—0.9 and the value of W\ from —0.6 to 0.6, we repeatedly calculate £i at z = 0.5, setting the 
cutoff mass scale at M c = 3.35 x 10 13 h~ 1 M & (the lowest cluster mass in the MICE catalog). 
The contour curves of £i in the Wq-w i plane are plotted in Figure [6l When wo is fixed, the 
relative abundance of the isolated clusters decreases as the value of W\ increases, which is 
consistent with Figure [51 

We also examine the degeneracy between Q m and W\ and between <7 8 and w% at z — 0.5, 
which are plotted in the left and right panels of Figure El respectively. For these plots the 
value of Wq is set at —1. As can be seen, when w\ is fixed, the relative abundance of the 
isolated clusters, £i, increases as the value of Q m increases. This implies that if there is 
stronger gravitational effect due to the larger amount of dark matter, the formation of a 
cluster can occur more easily even in isolated low-density environments. For a given value 
of £i, W\ increases with Q m . A similar trend is found in the ag-Wi degeneracy curves. When 
Wi is fixed, £i increases as the value of dg increases. The overall high density powers make 
it less hard for a cluster to form in the isolated low-density environments. 

It is worth mentioning here that the degeneracy trends sho wn in Figures EH3 are dif- 



ferent from those obtained from the total cluster abundance (see lWang fc Steinhardtlll998 



Weller et al.l 120021 ) . The increase of w\ has the same effect on the cluster abundance as the 
increase of Q m and <r 8 : the cluster abundance increases as the three parameters increase. 
In contrast, when the relative abundance of isolated clusters is used, the increase of w\ has 
the opposite effect: The relative abundance of the isolated clusters increases as Wi decreases 
and as fl m and cig increase. This result implies that the relative abundance of the isolated 
clusters may be helpful to break the degeneracy between W\ and the other key cosmological 
parameters. 



5. SUMMARY AND DISCUSSION 

We have derived the mass function of the isolated clusters in the frame of the recently 
developed Corasaniti-Achitouv formalism, assuming that for the case of the isolated clusters 
there is no disturbance from the surroundings and no ambiguity in the mass determination. 
The numerical results from the MICE simulations have been used to determine empirically 
the value of the drifting average coefficient, /3, which quantifies the degree of the deviation 
of the collapse barrier for the formation of the isolated clusters from the standard average 
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value. 

Extrapolating the validity of the Corasaniti-Achitouv formalism to the dynamic dark 
energy models and using our analytic result on the mass function of the isolated clusters, we 
have shown that the relative abundance of the isolated clusters, £j, defined as the ratio of 
the cumulative mass function of the isolated clusters to that of the non-isolated clusters at 
a given epoch, depends sensitively on the dark energy equation of state. This result proves 
our theoretical concept that the relative abundance of the isolated clusters is in principle a 
powerful probe of dark energy 

Yet, several additional works have to be done before constraining the dark energy equa- 
tion of state with £j. First, it will be necessary to derive the functional form of the drifting 
average coefficient, f3(z), for the isolated and non-isolated cases, separately. Second, it has 
to be confirmed if the Corasaniti-Achitouv formalism indeed works not only in a ACDM 
cosmology but also in dynamic dark energy models and how the drifting average coeffi- 
cient changes with the background cosmology. Third, it has to be examined whether or not 
the higher-order perturbation terms in the Corasaniti-Achitouv formalism have any non- 
negligible effect on the relative abundance of the isolated clusters. We plan to conduct these 
works and wish to report the results elsewhere in the future. 

We would also like to mention that our analytic model for the mass function of the 
isolated clusters will be useful not only as a probe of dark energy equation of state but also 
for predicting more accurately the cluster-related statistics since the isolated massive clusters 
may differ from the non-isolated ones in their physical properties such as average shape, 
dynamical state, gas entropy profile and etc (in private communication with E.Komatsu). 

I thank a referee for helpful suggestions and E. Komatsu for stimulating discussion. 
I acknowledge the use of data from the MICE simulations that are publicly available at 
:/ /www.ice.cat/mice| I also acknowledge the financial support from the National Re- 
search Foundation of Korea (NRF) grant funded by the Korea government (MEST, No. 2011- 
0007819) and from the National Research Foundation of Korea to the Center for Galaxy 
Evolution Research. 
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Fig. 1. — Mass function of cluster halos (top panel) and the ratio of the analytic prediction 
to the numerical result (bottom panel) at z — 0, 0.5 and 1 (in the left, middle and right 
panels, respectively). In each of the top panels the square dots represent the numerical result 
from the MICE simulations while the solid line is the analytic prediction based on the CA 
formalism. The errors are the one sigma scatter among the eight Jackknife resamples. 
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Fig. 2. — Best-fit value of the diffusion coefficient, Db, versus the linkage parameter, 6, used 
in the FoF algorithm to classify the clusters into the isolated and the non-isolated one at 
z = and 0.5 as dashed and solid lines, respectively. The result at z — 1 is very similar to 
the result at z = 0.5. 
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Fig. 3. — Same as Figured] but for the isolated cluster halos. 



- 18 - 




Fig. 4. — Ratio of the mass function of the isolated clusters to that of the non-isolated 
clusters at z — 0, 0.5 and 1 in the left, middle and right panel, respectively. 




Fig. 5. — Relative abundance of the isolated clusters (Eq. [15] ) for five different cases of the 
dark energy equation of state: w(z) = wq + W\zj (1 + z). The errors for the ACDM-case are 
obtained as the one a scatter among the eight Jackknife resamples. 
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Fig. 6. — Contour plot for the relative abundance of the isolated clusters at z — 0.5 in the 
Wq-Wi plane. 




Fig. 7. — Contour plot for the relative abundance of the isolated clusters at z — 0.5 in the 
Q m -Wi plane (left panel) and in the a 8 -Wi plane (right panel). 
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Table 1. Redshift, total number of clusters, and best-fit values of the two coefficients in 

the CAll formalism. 



z N T D B (3 



28i903i 0.38 0.107 
0.5 1684018 0.38 0.095 

1 749614 0.38 0.078 
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Table 2. Redshift, number of isolated clusters, and best-fit values of the two coefficients in 

the CA formalism 



z iV : D B (3 



1334200 0.0 0.34 
0.5 782254 0.0 0.20 

1 358073 0.0 0.08 



